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Discovery by Cassini's plasma instrument of heavy positive and negative ions within Titan's upper 
atmosphere and ionosphere has advanced our understanding of ion neutral chemistry within Titan's 
upper atmosphere, primarily composed of molecular nitrogen, with ~ 2.5 % methane. The external 
energy flux transforms Titan's upper atmosphere and ionosphere into a medium rich in complex 
hydrocarbons, nitriles and haze particles extending from the surface to 1200 km altitudes. The energy 
sources are solar UV, solar X-rays, Saturn's magnetospheric ions and electrons, solar wind and shocked 
magnetosheath ions and electrons, galactic cosmic rays (GCR) and the ablation of incident meteoritic 
dust from Enceladus’ E-ring and interplanetary medium. Here it is proposed that the heavy atmospheric 
ions detected in situ by Cassini for heights >950 km. are the likely seed particles for aerosols detected 
by the Huygens probe for altitudes <100 km. These seed particles may be in the form of polycyclic 
aromatic hydrocarbons (PAH) containing both carbon and hydrogen atoms C„H. V . There could also be 
hollow shells of carbon atoms, such as C B o, called fullerenes which contain no hydrogen. The fullerenes 
may compose a significant fraction of the seed particles with PAHs contributing the rest. As shown by 
Cassini, the upper atmosphere is bombarded by magnetospheric plasma composed of protons, HJ and 
water group ions. The latter provide keV oxygen, hydroxyl and water ions to Titan’s upper atmosphere 
and can become trapped within the fullerene molecules and ions. Pickup keV N-j, N* and CHJ can also be 
implanted inside of fullerenes. Attachment of oxygen ions to PAH molecules is uncertain, but following 
thermalization 0+ can interact with abundant CH 4 contributing to the CO and C0 2 observed in Titan’s 
atmosphere. If an exogenic keV 0+ ion is implanted into the haze particles, it could become free oxygen 
within those aerosols that eventually fall onto Titan’s surface. The process of freeing oxygen within 
aerosols could be driven by cosmic ray interactions with aerosols at all heights. This process could drive 
pre-biotic chemistry within the descending aerosols. Cosmic ray interactions with grains at the surface, 
including water frost depositing on grains from cryovolcanism, would further add to abundance 
of trapped free oxygen. Pre-biotic chemistry could arise within surface microcosms of the composite 
organic-ice grains, in part driven by free oxygen in the presence of organics and any heat sources, 
thereby raising the astrobiological potential for microscopic equivalents of Darwin’s “warm ponds" 
on Titan. 

Published by Elsevier Ltd. 


1. Introduction 

The Voyager 1 Titan flyby on November 12, 1980 showed a very 
complex encounter with Saturn’s magnetosphere with an induced 
magnetotail (Ness et al„ 1981) and a bite-out in the magneto- 
spheric electron population for £>500eV (Bridge et al„ 1981; 
Hattie et al„ 1982). It was also shown that the upstream plasma 
was composed of light and heavy ions (Hartle et aL, 1982; Sittler 
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et aL, 2004, 2005) and that the heavy ions (N + or 0 + ) had keV 
energies. The upstream magnetic field of -5 nT resulted in a heavy 
ion gyro-radii of -5000 km, which is larger than the size of Titan. 
A review of Cassini plasma and field observations of Titan’s 
interaction is given by Coates (2009a). The atmosphere detected 
by Voyager was primarily composed of N 2 , (Broadfoot et al„ 198 1 } 
as predicted by Hunten (1977) and the small methane component 
was detected earlier by Kuiper (1944). Voyager also detected a 
complex suite of hydrocarbons and nitriles using the infrared 
spectrometer IRIS (Hanel et al„ 1981, 1982). From these measure- 
ments exospheric models with H, H 2 , N, CH 4 and N 2 were 
constructed (Hartle et a!., 1982; Sittler et a!., 2004, 2005) and 
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complex neutral and ion chemistry was predicted to occur in its 
upper atmosphere and ionosphere (Yung et al., 1984; Yung, 1987; 
Toublanc et al„ 1995; Cravens et al., 1997). Sittler et ai. (2004, 

2005) also predicted that the dominant pickup ion slowing the 
magnetospheric plasma flow is CH 4 . The initial Cassini T A flyby 
revealed a similar interaction to that of VI but with the upstream 
heavy ion component uniquely identified as water group ions all 
of which contain oxygen. They also confirmed that the dominant 
pickup ion was CH 4 with some N 2 (Hartle et al., 2006a,b) as 
predicted (Sittler et al., 2004, 2005). The ion neutral mass 
spectrometer (INMS) instrument on Cassini during T A encounter 
also detected a very complex neutral composition in Titan’s upper 
atmosphere (Waite et al., 2005). During the T5 flyby INMS also 
detected a complex set of hydrocarbons which implied there was 
considerable ion chemistry occurring in this region (Cravens et al„ 

2006) . The rich ion hydrocarbon chemistry extended up to ion 
mass -100 amu. For later flybys when the spacecraft altitude 
dropped to 950 km Waite et al. (2007) reported the presence of 
benzene molecules plus heavy positive and negative ions with 
even higher masses. Coates et al. (2007) discovered heavy 
negative ions with mass -4000 amu and higher. The number of 
heavy negative ions increased with decreasing altitude down to 
950 km, the minimum altitude sampled by the Cassini spacecraft. 
Waite et al. (2007) and Coates et al. (2007) raised the possibility 
that the heavy ions were polycyclic aromatic hydrocarbons 
or PAHs. 

The above discoveries pose the question why are we seeing 
heavy negative ions or PAHs at all, In the studies of interstellar 
clouds the presence of such molecules has a long history. Herbst 
(1981) first posed the question that negative molecular ions could 
exist in molecular clouds. In that paper the high electron affinities 
(>2eV) for neutral free radicals and the large molecular size led 
to negative ion formation. Radiative attachment rate coefficients 
at cloud temperatures -10-50 K can be as high as the collision 
rate fc— 10~ 7 cm 3 /s. Later, Petrie and Bohme (2007) pointed out 
that there was no observational evidence for negative ions within 
molecular clouds, though their presence would be expected at low 
abundances. Subsequently McCarthy et al. (2006) combined radio 
observations of the star 1RC+10 and molecular cloud TMC-1 and 
their recent laboratory measurements to detect the negative ion 
C 6 H". This was later confirmed based on chemical models (Miller 
et al., 2007), which was then followed by detection of C S H~ and 
C 4 H - in the molecular cloud TMC-1 (Briinken et al., 2007). C 8 H~ 
was detected within the circumstellar envelope 1RC+10 216 by 
Remijan et al. (2007). Sakai et al. (2007) detected C 6 H within a 
low-mass star-forming region of L1527 and then later Sakai et al. 
(2008) reported the detection of C 4 H~ near the low-mass Class 0 
protostar IRAS 04368+2557, The anion C 3 N~ was discovered 
within the interstellar medium (ISM) by Thaddeus et al. (2008). 
Then Herbst and Osamura (2008) performed quantum mechanical 
model calculations confirming on theoretical grounds the likely 
presence of the detected Ion C„H~. They calculated high radiative 
attachment coefficients > 5 x 10~ 8 cm 3 /s that were similar to 
those measured by Miller et al. (2007). Relevant to the observa- 
tions reported here, the observed negative ions C„H~ have high 
carbon to hydrogen ratios (n>3 with n = 3 having the lowest rate 
coefficient and n = 8 the highest rate coefficient (Herbst and 
Osamura, 2008). 

In the case of Titan the ethynyl radical (C 2 H) and benzene 
radical C 6 H S (phenyl) can be important for the polymerization of 
PAHs (Goulay and Leone, 2006). At Titan the Cassini INMS 
instrument detected both acetylene (Waite et al., 2005) and 
benzene (Waite et al„ 2007; Vuitton et al, 2009a) and their 
radicals C 2 H (i.e., indirectly detected by Vuitton et ai., 2009b) and 
C S H S or phenyl by Vuitton et at. (2009a). Vuitton et al. (2007) also 
reported the detection of numerous hydrocarbons such as 


polyynes (C 4 H 2 , CiiH 2 , C g H 2 ), possible detection of methylpolyynes 
(CH 3 C 4 H, CH 3 C 6 H) and benzene (C S H 6 ). They also detected 
nitrogen bearing molecules cyanopolyynes (HC 3 N, HC S N), NH 3 , 
methanimine (CH 2 NH), nitriles (C 2 H 3 CN, C 2 H 5 CN) and possibly 
methylcyanopolyynes (CH 3 C 3 N, CH 3 C S N). Nitrogen bearing mole- 
cules can contribute to PAH formation and could thus be 
an important constituent of aerosols. The reader is referred to 
a comprehensive review article on PAH polymerization by 
Allamandola et al. (1989) the details of which are beyond the 
scope of this paper. The combination of recent detections of 
negative ions within molecular clouds and recent detections by 
Cassini of heavy negative ions at Titan suggests that there might 
be a connection between heavy negative ions and PAHs and/or 
fuilerene formations. This does not mean that neutrals and 
positive ions do not also play an important role in the growth of 
heavier molecules since reaction rates for ion-neutral collision, 
involving either positive or negative ions, tend to be much faster 
than those for neutral- neutral collisions (e.g., Waite et al., 2007). 
However, as noted above, the radicals C 2 H and C 6 H 5 are probably 
also present at Titan, in which case growth of fullerenes and PAHs 
can also occur via neutral-neutral collisions (Goulay and Leone, 
2006; Mebel et al., 2008; Bettens and Herbst, 1995). It is also true 
that acetylene and benzene can both lead to the growth of 
fullerenes and PAHs depending upon the chemical pathway taken 
(see Bettens and Herbst, 1995), The introduction of keV oxygen 
ions from Saturn's magnetosphere and their possible incorpora- 
tion into fullerenes and possibly PAHs as discussed here, adds 
another dimension with strong exobiological implications. The 
present paper explores the implications of these two potentially 
related discoveries, upper atmospheric formation of the heavy 
negative ions and bombardment by keV oxygen ions for chemistry 
of Titan’s atmosphere and surface. 

2. Observations 

2.1 Observations of heavy ions by Cassini 

We begin this section by reviewing the relevant plasma 
observations of Titan's upper atmosphere and ionosphere made 
by the Cassini plasma spectrometer (CAPS) and ion neutral mass 
spectrometer (INMS) instruments. The first discovery of extremely 
heavy negative ions within Titan's ionosphere was made by the 
high time resolution observations of the CAPS electron spectro- 
meter (ELS) during the T A encounter and revealed a surprising 
population of heavy negative ions. These ions were subsequently 
observed on 15 other Titan encounters, and were analyzed by 
Coates et al. (2007) and initially reported by Waite et al. (2007). In 
all cases the negative ions observed near closest approach, were 
narrowly confined to the spacecraft ram direction, and contained 
distinct energy peaks (see Figs. 1 and 2 from Coates et al„ 2007). 
During the spacecraft's flight through Titan’s relatively cold 
ionosphere the spacecraft motion, -6 km/s, provides an effective 
mass spectrometer feature for cold ionospheric ions. The conver- 
sion factor from ram energy to mass for singly charged ions at 
6.0 km/s is m amu = 5.32E CT (Coates et at., 2007). 

The CAPS instrument also has an ion beam spectrometer (IBS) 
primarily designed to measure high Mach plasma flows in the 
solar wind. But this instrument also provides good measurements 
of the positive ion component of Titan’s ionosphere which, during 
these encounters, also appear as a high Mach flow in the 
spacecraft reference frame. This instrument has detected heavy 
positive ions (Waite et al., 2007). The mass of the positive ions 
extends up to several hundred amu but the detected ions are 
significantly less massive than the detected negative ions. 
Abundances of both the positive and negative ions peak at the 
lowest altitudes encountered (-950 km). 
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For T a the observed maximum energy of the negative ions was 
~60eV, corresponding to an ion mass/charge m/q~320 amu/q. 
During other encounters, such as T16, the ion masses were as high 
as 10,000 amu/q as seen in Fig. 1 adpated from Coates et al. (2007). 
On various other encounters the ions were observed in rough 
mass groups at 10-30, 30-50, 50-80, 80-110, 110-200, (200-500, 
500+) amu/q (Coates et al„ 2007) because neither the ELS or IBS 
measure mass, but rather mass/charge. It was argued that these 
ions may in fact be multiply charged in this region of high electron 
density; therefore the actual mass could be up co five times 
higher. Such large highly charged ions provide seeds for forming 
aerosols (Coates et al., 2007, Waite et al, 2007), This relationship 
between mass and charge state is strongly tied to the assumed 
size of the ions with the following relationships. Waite et al. 
(2007) used cp — 2.5kT/e -0.027 V (r~125l<) based on the 
surface potential for a small dust grain with its corresponding 
charge Q = 4m 0 a(i>exp(-alA D ) (e 0 is permittivity of free space, a is 
grain's radius and /. D is the Debye length). As discussed below, if 
the heavy ions are fuilerenes, then the mean ion radius is 
estimated to be r A ~12A giving Q~1 and a maximum mass 
< 10,000 amu. If one uses for the heavy ions PAHs, then r A ~60 A is 
possible so that Q~5 is estimated and very large ion masses 
become possible. Therefore, whether the heavy negative ions are 
fuilerenes or PAHs becomes an important discriminator with 
regard to their true size. 

As discussed in Coates et al. (2009b), negative ions were not 
anticipated this high in Titan's atmosphere and were not included 
in pre-Cassini chemical schemes (e.g., Wilson and Atreya, 2004) so 
that new chemical models are required. As noted earlier, ion- 
neutral collisions have high collision cross-sections, which in the 
case of the heavy ions approach the geometric cross-section 
(Waite et al., 2007). This allows such ions to grow rapidly in mass 
and size resulting in heavy ion formation. An important parameter 
in this growth process is the charge state. In case of negative ions 


this depends on the electron affinity of the molecule which as 
discussed below favors fuilerenes. 

The relationship between these ions and the heavy positive ion 
population was discussed by Waite et al. (2007, 2008). They 
suggested that nitrogen and methane in Titan's high atmosphere 
would be acted on by sunlight and magnetospheric particles 
forming heavier but relatively simple species by dissociation and 
ionization processes. Eventually this could lead to the growth of 
benzene and the other heavier ions detected by IBS up to 
—350 amu/q (note that the upper limit could be set by the IBS 
detection threshold and even heavier positive ions could exist in 
substantial amounts). The growth process would continue to form 
heavier positive and negative ions. They suggested that these may 
be the tholins originally postulated by Sagan and Khare (1979), 
With continued growth these molecules could become large 
enough to be called aerosols and drift down towards Titan's 
surface. This chain of processes, initiated by magnetospheric 
interactions at the top of the atmosphere eventually affects the 
surface composition of Titan. This idea was supported by the 
extensive observations of negative ions by Coates et al. (2007). 

Recently, further analysis of negative ion signatures has 
extended the number of encounters where negative ions were 
observed to 22 during Cassini's prime mission (Coates et al., 
2009b). This allows a systematic study showing that the higher 
mass negative ions are observed preferentially as follows: (1) at 
low altitudes, with the highest mass (~10,000 amu/q) at Cassini’s 
lowest altitude of 950 km, (2 ) at high Titan latitudes and (3 ) in the 
region of the terminator. On the basis of these results, Coates et al. 
(2009b) suggested that the formation of high mass negative ions 
is more efficient, or that destruction is less efficient when solar 
flux is highly attenuated. Coates et al. (2007) also suggested the 
possible ion identifications: (1) 10-20 amu/q CN , NH 2 and 0~; 
(2) 30-50 amu/q NCN~, HNCN and C :J H-; (3) 50-80amu/q C s Hg, 
CgHT C fi Hi ; and (4) mass groups 80-110, 110-200, 200-500 and 


Energy/charge (eV/q) 

1 10 100 1000 



Flg.l. PlotofCAPS-ELS ion counts versus ion energy (keV) or ion mass-per-charge (M/Q) from Coates et ai. (2007; rop panel and same for bottom panel except here we have 
ion density. 
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500-10,000 as polyynes, high order nitriles, PAHs and cyan 
aromatics. Nitrogen substitution PAHs are another possibility, but 
may not be stable in the UV environment of Titan’s upper 
atmosphere. They may, however be stable in the haze layer of 
Titan's stratosphere at altitudes below 400 km (e.g., Ricca et al„ 
2001 ). More recently, Vuitton et al. (2009b) using a photochemical 
model have identified peaks in the CAPS-ELS mass spectrum at 
22±4amu/q as CN~/C 2 H - , 44±8amu/q as C 3 N“/C 4 H~ and 
82-14amu/q as C 5 N~. Most noteworthy is C 2 H~ which is more 
abundant than CN~ below 850 km altitude in their model 
calculations. The fact that the ionization ratio for C 2 H~/C 2 H is 
~2.4 x 10 5 in their model shows that the free radical C 2 H is likely 
present in Titan's upper atmosphere. 


2.2. The formation of fullerenes and PAHs in Titan's upper 
atmosphere 

As discussed in the papers by Coates et al. (2007, 2009b), and 
Waite et al. (2007, 2008) the heavy positive and negative ions 
could be evidence for the formation of PAHs (C 2n H x ) in Titan’s 
ionosphere. The heavier ions are observed at lower altitudes with 
9501<m being the lowest altitude sampled by the Cassini space- 
craft. Vuitton et al. (2009a) have presented evidence of benzene in 
Titan’s upper atmosphere with molar fractions ~10~ H at 950 km 
altitude and levels higher than expected (see Lawas et al., 
2008a, b). Acetylene C 2 H 2 has a molar abundance ~10~ 3 at 
1100 km (Vuitton et al„ 2009b) and is more abundant than 
benzene (also see Shemansky et al., 2005; Waite et al„ 2007). For 
atmospheric temperature like that at Titan, 7M50K, one expects 
acetylene and benzene to lead to the formation of fullerenes and 
PAHs, respectively. 

Laboratory measurements at much higher temperatures can 
produce fullerenes using low pressure benzene-oxygen flames 
(Gerhardt et al„ 1987; Richter and Howard, 2000). PAHs can be 
formed in the laboratory using low pressure acetylene-oxygen 
flames (see Gerhardt et al., 1987). As discussed below at lower 
temperatures T<300K the reverse is true. There is a lot of free 
energy within Titan’s upper atmosphere from bombarding high- 
energy photons (UV/EUV and X-rays) and magnetospheric 
charged particles. This high-energy input results in high thermal 
electron temperatures r e< ~1000K (Wahlund et al., 2005) and hot 
suprathermal photoelectrons T e h~20,000 K (Hartle et al., 2006a,b; 
Coates, 2009a). But, the ion-neutral and neutral-neutral collisions 
that lead to the growth of the heavy ions occur for T~150K. The 
high electron temperatures may affect ion charge state and 
formation of radicals. The above measurements were also made in 
oxygen flames, but Titan is essentially devoid of oxygen except for 
low abundances of CO, C0 2 and H 2 0 within Titan's atmosphere. 
Therefore, we feel these earlier laboratory measurements are not 
applicable to Titan. 

More relevant to the observations at Titan fullerene and PAH 
formation occur under space environmental conditions within 
interstellar and circumstellar environments that more closely 
resemble Titan’s upper atmosphere (e.g„ Bohme, 1992). X-ray 
crystallography measurements have shown fullerenes to be 
hollow carbon cages (see Fig, 2) made of 60 carbon atoms C 60 
(720 amu) which have spherical shape or 70 carbon atoms C 70 
(840 amu) with ellipsoidal shape (Kratschmer, et al., 1990). The 
measured inner diameters are ~7 A with effective mass density 
pi.-~l.8gm/cm 3 (i.e., mean radius r~5.4A for p~2.25 gm/cm 3 for 
outer shell made of carbon with thickness ~3.4A). This should be 
compared to aerosol radii r A ~12.5nm estimated by Liang et a!. 
(2007). Note, the r A --260nm estimated by Waite et al. (2007) was 
an over estimate (Waite, 2008, private communication). Both C«o 
and C-o are known to be very stable (Kroto et al., 1985; Kroto et al., 
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Fig. 2. Shows cage structure forCeo fullerenes as described in the review article by 
Bohme (1992;. (Three examples of proposed "winged" structures for double 
adducts of the type CsofXH)!' where XH=NH 3 ,(CH 3 ) 2 CO,and CH 3 CN. 


1991) and are found to be the end product of forest fires (Richter 
and Howard, 2000; Calcote et al., 1990; Frencklach, 2002). The 
main difference between fullerenes and PAHs is that the former 
has no hydrogen (C 6 o). while PAHs have hydrogen (C 2 , t H x ). 

There is a "zipper” mechanism often discussed with regard to 
the formation of fullerenes. Hunter et al. (1994) gave one 
explanation on how this “zipper” mechanism works. They based 
their conclusions on laboratory measurements where they used a 
laser vaporization of graphite source which produced ionized 
carbon chains and rings. The cluster ions of carbon were injected 
at different energies into a Helium buffer gas where they were 
thermalized. They use a quadrupole mass spectrometer to select 
cluster ion products of different mass. Their measurements 
showed that fullerenes were produced from intermediate mono- 
cyclic, bicyclic and tricyclic (i.e., polycyclic) carbon rings. Then 
from low frequency vibration modes and strain relief of the larger 
cyclic rings and the resulting spiraling of polyyne chains, 
fullerenes form. These measurements also showed that the 
activation energy ~2.4eV for fullerene formation were unexpect- 
edly low. In the papers by Thaddeus (1994, 1995) the argument is 
made that PAHs are not likely to be the dominant path to heavy 
molecules or aerosols and grains within diffuse molecular clouds, 
but rather fullerenes which have their origins from long carbon 
chains which then fold into polycyclic rings and then fullerenes. 
Thaddeus (1994, 1995) then go on to say that fullerenes from C + 
insertion reactions with H, 0, N, S, etc, can then grow into larger 
“amorphous carbon" organic polymers. There is experimental 
evidence of clustering of fullerenes (C 6c with n < 147 observed in 
the laboratory by Martin et al. (1996) with major peaks in their 
mass spectrum at n = 13 and 55. Whether such clustering within 
Titan’s atmosphere can occur must await further laboratory 
measurements more characteristic of Titan’s atmosphere. 
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If heavy negative ions observed by Cassini are fuilerenes, then 
they should prefer a negative charge. Laboratory measurements 
by Yang et al. (1987), Wang et al. (1991) and Wang et al. (2005) 
have shown C 60 to have high electron affinity (2.6-2.9eV), and 
since electron attachment for PAHs occur for electron affinities 
>leV (Moustefaoui et al., 1998), fullerene negative ions are 
allowed. But, in some cases fuilerenes can be stripped of electrons 
to form +1 and +2 cations (Weiske et al„ 1991 ). The measurements 
by Moustefaoui et al. (1998) were made at relatively high electron 
temperatures ~300K or more and that it is difficult to achieve 
lower electron temperatures T~-150K in the laboratory (Coulay 
et al„ 2004). Titan's upper atmosphere and ionosphere are 
relatively cold with neutral and ion temperatures -150K (Waite 
et al., 2005, 2007), but the electron temperatures can be quite hot 
with thermal electron temperatures T e , c ~1000K (Wahlund et al., 
2005) and hot component photoelectron temperatures 
T e- i,~20,000 K (Hartle et al., 2006a, b; Coates, 2009a; Cravens 
et al., 2009). It is important to note that electron attachment to 
fuilerenes C GQ and C 70 has an exponential temperature depen- 
dence that increases rapidly for T e > 500 K with relatively high rate 
coefficient /? e ~3 x 10~ s cm 3 /s (Smith et al., 1993; Spanel and 
Smith, 1994; Spanal and Smith, 1995; Jaffke et al., 1994) at 
T ec ~1000K (Wahlund et al., 2005), so electron attachment to 
fuilerenes seems allowed. More recently, Agren et al. (2009) did a 
more systematic study of Titan's ionosphere and determined the 
dayside peak densities to be 2500-3500 el/cm 3 and the nightside 
peak densities to be 400-1000 el/cm 3 . The ionospheric peak 
electron temperature is rather constant near 350-700 K. The 
lower electron temperatures given in this report will give lower 
electron attachment rates for C so ~10 s cnt 3 /s and for electron 
densities N ec ~3000 el/cm 3 on the dayside attachment time 
scales are ~9 h, while for nightside with N ec ~ 1000 el/cm 3 electron 
attachment time scales ~27h. If one uses the temperature 
for photoelectrons with r e h~20,0001< and N eh -~10 el/cm 3 (Hartle 
et al.. 2006b for TA) the fullerene attachment rate is estimated 
to be /f e h~4.4 x 10~ 7 cm 3 /s with attachment time scales r eh ~l/ 
(N, ? h/?c.h} ~2.6 days which are longer than the thermal electron 
attachment time scales. 

In case of positive ions recombination time scales are a~10 -5 / 
Tfi 6 ~2 x 10~ 7 cm 3 /s for T CC ~700K are relatively fast T ae '~l-2h 
(Cravens et al., 2006). Photoionization rates for atoms and smaller 
molecules, for which the ionization potentials are >10eV, are 
typically ~10 s s (see Sittler et al., 2008) or time scales ~3yr. 
Borucki et al. (2006) considered photoemission rates for embryo 
aerosols (3 A) within Titan's lower atmosphere for heights less 
than 400 km. Fullerene sizes are similar to the embryo aerosols 
considered by Borucki et al. (2006) and have ionization potentials 
-7.07 eV (Bettens and Herbst, 1995), so photoemission rates for 
fuilerenes and their clusters are expected to be higher than those 
for atomic and molecular ions. Their calculations do show that 
when in shadow photo-detachment of electrons from negative 
ions (i.e., fuilerenes) will be much lower and thus more likely to 
retain their negative charge. This is supported by the CAPS-ELS 
observations reported by Coates et al. (2009b) where heavy 
negative ions are more important near terminator local times at 
Titan. Estimation of the actual charge state of fuilerenes requires 
development of a complex model similar to that by Borucki et al. 
(2006), including photo-detachment rates, ion-neutral and ion- 
ion charge transfer rates and is thus beyond the scope of this 
paper. But, the relatively high electron affinities and electron 
attachment rates are consistent with negative charge states for 
fuilerenes. This will be the subject of future work. 

Ionization of the upper atmosphere comes in part from solar 
UV, magnetospheric electrons and X-ray photons but also arises 
from bombardment by hot keV magnetospheric heavy oxygen 
ions (Hartle et al., 2006a,b,c) and energetic protons and oxygen 


ions (Mitchell et al., 2005 and Cravens et al., 2008). The energetic 
protons and >10keV oxygen ions can penetrate below 9501<m 
(Luna et al., 2003; Cravens et al., 2008; Hartle et al„ 2006c) where 
the heavy positive and negative ions are observed by Cassini 
(Coates et al., 2007; Waite et al., 2007). Galactic cosmic ray 
protons (Capone et al., 1983; Molina-Cuberos et al., 1999) 
penetrate much more deeply and dominate ionization in the 
lower atmosphere. Ionization layers at altitudes -700 km can 
form due to meteoritic dust impacts as shown by Molina-Cuberos 
et ai. (2001 ). Ionization layers -500 km were observed by the 
Cassini radio science experiment by Kliore et al. (2008). The key to 
the formation of fuilerenes and PAHs within this environment is 
then the combination of relatively high methane abundance and 
the high-energy radiation input to the atmosphere. 

2.3. Role of oxygen input from magnetosphere and micrometeorites 

The role of oxygen becomes particularly significant for 
astrobiology if fullerene formation occurs. Oxygen is an important 
ingredient of amino acids and other chemical compounds 
necessary for evolution of life, and exothermic oxidation of 
nutrients is critical to biochemical processes. Major abundance 
of free oxygen in Titan's highly reducing atmosphere is highly 
unlikely (Raulin and Owen, 2003), but the bound molecular form 
is found in trace species such as CO at the 60 ppm level in Titan’s 
atmosphere as discovered by Lutz et al. (1983). This was later 
followed by the discovery of CO 2 in Titan's atmosphere at the 
1.5 ppb level by Samuelson et al. (1983), who showed that water 
entry from above may be required in order to explain the presence 
C0 2 having a relatively short lifetime -5 x 10 4 yr. In both cases the 
oxygen is locked up in the CO and C0 2 which are fairly stable 
molecules ( CO - C+O+l 1.17 eV; C0 2 - CO+O+5.46 eV; C0 2 - C+O+O+ 
16.63 eV). It was argued that the CO was probably primordial in 
origin but the recent model results by Horst et al. (2008) showed 
that this was not required due to the entry of keV oxygen ions to 
Titan's upper atmosphere (Hartle et al., 2006a, b). 

Oxygen can be introduced into Titan’s upper atmosphere 
by entry and ablation of micrometeorites containing water 
molecules (Samuelson et al., 1983). Models of Titan’s atmosphere 
after the discoveries of CO and C0 2 required an inward flux 
of oxygen as summarized in Table 1. In the case of Yung et al. 
(1984) they required oxygen influx -6.1 x!0 s O atoms/cm 2 /s, 
while that by Toublanc et al. (1995) used a downward flux 
—1.5 x 10 6 O atoms/cm 2 /s. The mechanism of meteoroid ablation 
in Titan's upper atmosphere was initially suggested by ip (1990) 
for the case of water ice particles and stony particles. Then English 
et al. (1996) proposed a micrometeorite ablation model for water 


Table 1 

Summary of oxygen, source for Titan's upper atmosphere. 



Oxygen source 

In flux rate 
(0/cm 2 /s) 

Lutz et al. (1983) 

CO discovery -60 ppm 


Samuelson et al. 

C0 2 discovery -1.5 ppb ®» 1 10 mbars E 

6.1 >: 10 s 

(1983) 

ring & Meteor icand meteoric dust as 
source of O 


Yung etaf. (1984) 

Based on CO &at)d C0 2 chemistry 

6.1 X 10 s 

Toublanc et al. 

Chemical model with CO &and CDs* 

1.5 x 10 r - 

(1995) 

English et al. (1996) 

Micrometeorite dust ablation 

3 x 10‘ 

Lara et al. {1996} 

CO &arsd C0 2 model based on English 
et ai. (1996) 

3 x 10° 

Coustenis et al. 

Discovery of H 2 0 ~8 ppb @ 400 km 

0.8-2.8 x 10 e 

(1998) 

Haiti e et al (2006a) 

Ta keV magnetospheric 0‘ 

I.S x to'- 

Hartle et ai. (2006b) 

Ta keV magnetospheric O' 

1.1 x 10 s 
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entry with peak deposition rate at -700 km altitude. Their model 
showed peak flux -3 x 10 6 mol/cm 2 /s. The discovery of water in 
Titan's atmosphere by Coustenis et al. (1998) with mole fraction 
~8 ppb at 400 km altitude ended any doubts about the presence of 
water products in Titan's upper atmosphere. The Coustenis et al. 
(1998) results inferred an influx rate of 0.8-2. 8 x 10 6 mol/cm 2 /s. 

Cassini results continue to build the case for oxygen and water 
sources in the upper atmosphere. Hartle et al. (2006a, b) presented 
the first observational evidence of magnetospheric keV oxygen 
ions into Titan's upper atmosphere with an estimated influx rate 
~1.1 x 10 6 0 atoms/cm 2 /'s at the exobase. The presence of negative 
oxygen ions in Titan’s upper atmosphere can be inferred by the 
mass 16 peak in Fig. 1 from T16 adapted from Coates et al. (2007). 
This might be expected due to the high electro negativity of 
oxygen atoms. This is also suggested by Cravens et al. (2008). 
although Vuitton et al. (2009b) based on T40 heavy negative ion 
measurements argue against a significant presence of O”. When 
Vuitton et al. (2009b) included the influx of oxygen — 10 s 0/cm 2 /s, 
which is 5% of “typical" of that reported by Cravens et al. 
(2008) for T5 and 10% of that reported by Hartle et al. (2006a,b) 
~10 6 0/cm 2 /'s forTA, the O - densities increased from a maximum 
of 10 ' 6 Cr/cm 3 at 1100 km altitude to 4 x 10“ 3 0~/cm 3 at 1100 km. 
If a more typical value of the influx is used and the uncertainties 
in the cross-sections are accounted for one could argue that 
elimination of 0“ from the CAPS-ELS mass spectrum cannot yet be 
made. It is also important to note that there are times when no 
heavy ions are present in the upstream magnetospheric flow and 
times when it dominates the upstream flow (Sittler et al., 2009) 
and therefore one must consider such complications when 
interpreting negative ion mass spectra for specific Titan flybys 
(i.e., T16 vs. T40). 

The CDA observations of E-ring dust particles by Srama et al. 
(2006) within Saturn's outer magnetosphere allows one to 
estimate an influx rate of water molecules into Titan's upper 
atmosphere -2.5 x10 s water mol/cm 2 /s at the exobase. In the 
case of magnetospheric oxygen ions and E-ring dust particles, the 
source of the water and oxygen can be traced back to Enceladus 
(Johnson et al., 2005; Sittler et al., 2006). 


2.4. Trapping of free oxygen and hydroxyl ions in seed particles 

Laboratory experiments show that C|o 6~161<eV ions (lab 
frame kinetic energy) can trap lower mass target gas atoms inside 
the fullerenes. Magnetospheric keV 0 + ions incident upon Titan’s 
upper atmosphere can be implanted and become trapped within 
Titan's fullerenes. The laboratory measurements are quite ex- 
tensive for fullerene trapping of target noble gases such as He, Ne, 
Ar, Xe and !<r, Weiske et al. (1991) found that C 60 + He — C*He 
(46<x<60) for He target gas at Cg 0 incident energy (product 
fractional percentage) of 6keV (30%), 8keV (34%), lOkeV (14%) 
and 16 ke¥ (<1%). The equivalent center of mass energies for 
incident He ions would be 33, 44, 56 and 89 eV, respectively. For 
incident 0* ions the equivalent center of mass energy would range 
from 132 to 356 eV. As keV CT ions lose energy within Titan's 
upper atmosphere (e.g., Shah et al., 2009) and their energies drop 
below 400 eV they can become trapped in fullerenes with higher 
probability. This trapping will occur deeper in the atmosphere as 
the CT incident energy increases (i.e., ion penetration into 
atmosphere increases as the ion inititai energy increases). In case 
of Argon center of mass energies used were -450 eV, while for Kr 
center of mass energies were ~1230eV (Caldwell et al., 1991, 
1992). The higher center of mass energies can cause fragmenta- 
tion of the fullerene with ejection of carbon atoms. Higher center 
of mass collisions up to 1 keV for O’ have yet to be done, in some 
cases more than one He atom was implanted into a fullerene. 


Caldwell et al. (1991, 1992) did similar experiments using He, Ne 
and Ar and found end products C S0 He + , CeoHeJ, C 70 He + , C 8 4 He\ 
C ss Ne + and C 55 Ar + . Xe and Kr were not captured due to their larger 
size as compared to the space inside fullerenes. These results are 
consistent with X-ray crystallography measurements of fullerene 
dimensions and size within their hollow cage. 

From these laboratory results two mechanisms are invoked. 
Mechanism 1 is the “low pressure limit” for which the ion is 
injected directly into the fullerene, while mechanism 2 is in the 
“high pressure limit" for a collision induced fracture of the carbon 
cage occurring first and then being followed by atom capture 
as an ion complex. The first mechanism could apply at Titan 
for direct injection of oxygen ion into the fullerene. Although, 
most of the oxygen ions will form CO and C0 2 once thermalized, 
the direct implantation mechanism does not have this problem. 
The capture cross-section (fullerene geometric cross-section 
o-F~7rrF~4 x lCT ls cm 2 ) is quite large. Results of Coates et al. 
(2007) for mean negative ion mass — ■ 4000 amu (i.e., embryo 
aerosol) give 4000/720-6 cages, radius r A ~ 7 A or larger and 
geometric cross-section er A ~7iTA~1.5 x 10~ 14 cm 2 . 

Trapping of oxygen ions inside the fullerenes provides a stable 
means of oxygen conveyance within the upper atmosphere which 
is ultimately delivered to the surface. Complex chemistry of pre- 
biotic and even biotic interest could occur as oxygen-laden aerosol 
particulates precipitate to any surface hot spots associated with 
local cryovolcanism (Sotin et al, 2005: Lopes et al., 2007; Nelson 
et al., 2008) and where liquid water might be present near the 
surface. 

In the modeling by Michael and Johnson (2005) and Michael 
et al. (2005), the energy input into Titan’s upper atmosphere by 
keV magnetospheric heavy ions, the heavy ions tended to 
dominate the energy input to Titan's upper atmosphere. They 
assumed the ions were directed down into Titan's atmosphere. 
This was the case since for VI the convective electric field pointed 
radially away from Saturn so that pickup ions on the Saturn side of 
Titan would be accelerated into Titan’s upper atmosphere. For T5 
the electric field is pointed southward so pickup ions are 
accelerated into Titan’s northern hemisphere. In cases when Titan 
is above the magnetospheric current sheet the pickup ions will be 
accelerated into Titan's southern hemisphere. They found the 
pickup ions tended to dominate the energy input to Titan’s upper 
atmosphere relative to the magnetospheric input. Therefore, 
implantation of N + , N 2 and CHJ pickup ions into the fullerenes 
could also be occurring at levels comparable to oxygen input rates. 

2.5. Transport to lower atmosphere and surface as aerosols 

Estimates of the mass flux of the heavy ions observed in Titan's 
ionosphere downward to the surface can address whether these 
ions might account for the aerosols observed by Tomasko et ai. 
(2005) using the Huygens D1SR Imaging system. In principle we 
can also estimate the potential surface delivery rate of oxygen 
bound in fullerenes due to aerosol transport. As stated earlier the 
heavy negative ion density increased with depth in the atmo- 
sphere with the highest masses and number densities at the 
minimum altitude 950km probed by Cassini (Coates et al., 2007). 
Such increases likely continue to lower altitudes and perhaps 
ultimately to maximum densities of aerosol cloud layers 
100-200 km or lower (de Kok et al., 2007) above the surface. 
Lavvas et al. (2008a, b) developed a photochemical model of 
Titan's atmosphere which extended from the surface to the upper 
atmosphere z ~ 1400 km altitude. In these papers the size 
distribution of the haze particles was derived with height using 
a 1 D diffusive transport code with source and sink terms for the 
size distribution of particles. The calculations determined the size 
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distribution of the aerosols at different heights with peak number 
densities ~100 part/cm 3 and size r A -10~ 3 pm at 950 km. Lavvas 
et al. (2008b) give particle mass -5000 amu, very close to that 
observed by Coates et al. (2007) at 950 km. Near the surface these 
latter authors estimate particles sizes —0.5 pm with number 
densities —7 part/cm 3 . 

Borucki and Whitten (2008) and built on earlier models by 
Borucki et al. (1987, 2006), modeled haze particle distributions for 
heights less than 150 km. They assumed three mass fluxes 
4 x 10” l4 kg/m 2 /s for case B, 4 x 10" 13 kg/irri/s for case C, and 
10' 12 kg/m 2 /s for case D, also assume mass flux conservation, and 
do not use diffusion terms which are essentially unknown. Using 
Stokes Equation from Roe et al. (2002), they estimated fall velocity 
-5 x 10 _fl m/'s for case B and -4 x 10^ 5 m/s for case D. The particle 
sizes were r-0.45 pm for case B, 0.95 pm for case C and r~1.3 pm 
for case D. The particle abundances were —6 x 10 6 part/m 3 ~6 
part/cm 3 , which are very close to that estimated by Lawas et al. 
(2008b). In this paper we estimate a simplified fall velocity 
applicable to the higher altitudes probed by Cassini which 
includes simple frictional drag from molecular collisions on the 
heavy ions. The expression is as follows: 


Adrift — 


M se edg( Z) 
N N2 M Nj WN 2 <T seed 


(1) 


with 

ff seed j = nr seedj ® 


where M seca is the mass of the heavy ion observed by Coates et al. 
(2007) and we refer to as seed particles for the aerosols (or 
embryo aerosols), g(z) is the local acceleration of gravity which is 
a function of height z, N n , is the number density of N 2 molecules, 
M n , = 28 amu is the mass of N 2 , w Nj is the thermal speed of the 
N 2 molecules and cr seed is the geometric cross-section of the seed 
particles plus that of the N 2 molecules. We use this expression to 
estimate the heavy ion fall speed at 950 km which one finds to be 
-10 m/s for fuilerenes and -0.4 m/s for PAHs (i.e., both with same 
JVf/Q, see previous discussions about ion charge state and cross- 
sectional size). 


Negative Ion Fall Velocity 



Fig. 3. Plot of aerosol fali velocity (mis) versus altitude based on our simplified 
expression in Eqs. (1 ) and (2; from above 950 km down to the surface. Two curves 
are shown for particles mass 4000 and 40,000 amu (top curve). !f particles are 
spherical they will fall taster as they grow, but if more planar then their fell 
velocity is independent of particle mass for a given height. We also show the 
estimates by Borucki and Whitten (2008a,b) for altitudes < 150 km, Our estimate 
is less than Borucki arid Whitten ( 2008a, b) as would be expected, since the particle 
size ~1 ftm is >6.5 A (4000 amu; and 14 A (40,000 amu) for seed particles 
observed above 950 km altitude. Our expression for spherical 1 pm sized particles 
would yield fall velocities several orders of magnitudes greater than that shown 
here for 40,000 amu particles. 


We show in Fig. 3 the fall speed of the haze particles as a 
function of height for different assumed haze particles mass (4000 
and 40,000 amu). In reality, as the haze particles fall they will 
grow in size and tend to fall faster as they become more massive 
except for the increase in frictional drag due to the increased 
cross-sections and atmospheric N 2 density with depth. Looking at 
Eqs. ( 1 ) and (2), heavier spherical ions will fall faster linearly with 
particle radius r(i.e., M P ~(47t/3)pr 3 and a P ~jrr). Fluffy or planar 
particles will fall at a speed independent of particle mass at a 
given height. Assuming fuilerenes as the primary aerosol 
component, and taking the ion number densities with mean 
mass -4000 amu of the heavy negative ions we estimate at 
950km altitude a downward mass flux $—2.7 x I0~ 14 kg/m 2 /s. 
This flux is smaller than the lower limit $—4 x 10 _,4 kg/m 2 /s used 
by Borucki and Whitten (2008a,b) and constrained by Tomasko 
et al.’s (2005) aerosol observations but considering uncertainties 
identical. If the mass density flux increases below 950 km as 
expected, then the heavy ions observed above 950 km may be 
sufficient to account for all the aerosols observed by Tomasko 
et al. (2005) near the surface. If one uses the mass flux rate 
—4 x 10~ ,4 kg/m z /s and a surface mass density —420 kg/m 3 it will 
take — 300Myrs to accumulate an aerosol depth of one meter 
(420 kg/m 2 ) uniformly distributed around the globe. If one uses 
the higher mass fluxes of Borucki and Whitten (2008a, b) the 
accumulation times could be as short as — lOMyrs. 


2.6. Loss of ionospheric ions to Titan's surface 

We now consider the effect of heavy negative ion drift and 
formation of an ambipolar electric field. We have ignored the role 
of heavy positive ions N HeMy j since they are less abundant than 
the dominant ionospheric ions N I0 nj. Then assuming the seed 
particles are negative and imposing the charge neutrality 
condition 

, N S eaJjQseedj T Ng — ^ ] hflONJ T ^ ^ HeavyjQHeavyJ (3) 

; j 


Ne s £N,o N j - (4) 

j 

for which the heavy negative ions )C;NseedjCLeed,;<N e , impose zero 
current condition 

A/eed.;Qvvdjf(lrtft. i ,eed j ~ , Nj Qn0 j V)j r jftj onQ ; = 9 fa) 

J i 

(the electrons are assumed to be collisionally coupled to the 
neutral atmosphere) and frictional drag one can derive the 
following expression for the ambipolar electric field 

E(Z) — If Q>wd i /' AT-Z ) ~ /rtc)] g 

etNseedQLd + NlONfffseedM)] 

The following parameters are used (Tse«t/<Ts— 40 (typical elastic 
scattering cross-section is a s = 5 x 10~’ s cm 2 for ion-neutral colli- 
sions), N[ 0 n = 5000 ions/cm 3 (HCNH* ion), Ns «d = 200 anions/cm 3 , 
Msccd/MioN— 143 (Mseeci = 4000 amu), seed particle charge state 
Qseed = 1. E(z) = 0.180 pV/m, a^d = e£(z)/Mseed~0.43 cm/s 2 «$z) = 
72 cm/s 2 at 950 km, a [ON = e£(z)/M 10N -62 cm/s 2 — g(z) at 950 km, and 
then with Vd,, f ; S<H d— lOm/s and Vd riffi ono-0.4 m/s. Using these we 
estimate a loss of positive ionospheric ions to the lower atmosphere 
by negative ions dragging them along —2 x 10 3 ions/crrri/s (i.e., at 
lower altitudes they probably recombine and become part of the 
lower atmosphere). If, electrons were free to move vertically, they 
would tend to short circuit this electric field. The net effect would be 
to reduce die downward drag on the ionospheric ions by the heavier 
negative ions. This ioss is not very large, and about 4% of the 
ionospheric wind loss from the topside ionosphere (Hartie et ai„ 
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2008: Sittler et al., 2009). The atmospheric loss from the fall of 
negative ions in units of mass 28 is much larger: ~5 x 10 fa ions/cm 2 /s 
and when one combines the two losses the total ion loss ~5.2 x 10 B 
ions/cm 2 /s. This translates to a loss rate of 2.4 x 10 16 gm/cm 2 /s and 
can be compared to a loss rate ~5 x 10”' 6 gm/cm 2 /s at 1000 km 
altitude estimated independently by Wahlund et al. (2009). 

2.7. Other observations and final comments 
2.7.1. Fullerenes versus PAHs 

We first address the issue of fullerenes versus PAHs with the 
former composed only of carbon C,„ while PAHs are composed of 
both carbon and hydrogen C„H m . Hydrocarbons like C 2 H 2 are 
considered unsaturated in hydrogen, while C 2 H 6 is considered 
saturated. It is sometimes implied that since Titan’s atmosphere is 
rich in hydrocarbons with relatively large amounts of hydrogen, that 
hydrogen poor species such as fullerenes would not form. But, in 
fact the reverse may be true and many experts in the study of ISM 
such as Thaddeus (1994, 1995) favored fullerenes over PAHs in an 
environment similar to that for Titan (i.e., low temperatures 
T< 150 K). Observation of negative ions in the interstellar medium 
were primarily composed of polyyne radicals C 2n H~ (McCarthy 
et al., 2006; Brunken et al„ 2007; Sakai et al., 2007; Remijan et al„ 
2007) or unsaturated carbon nitrogen species such as C 3 N“ 
(Thaddeus et al., 2008). Vuitton et al. (2007) and later Vuitton 
et al. (2009b) showed significant abundances of poiyynes (C 2 H 2 , 
C 4 H 2 , C 6 H 2 , C s H 2 ) in Titan’s atmosphere at 1100 km using a 
photochemical model constrained by 1NMS ion measurements. In 
Vuitton et al. (2007) they also note that the abundance for the 
higher carbon chain poiyynes within Titan’s atmosphere did not 
decrease as fast as that in other photochemical models (Yung et al., 
1984; Toublanc et al., 1995; Wilson and Atreya, 2004), the ISM 
(MacLeod et al., 1984; Remijan et al., 2006; Cernicharo et al., 2001 ) 
and/or laboratory measurements (Coll et al., 1999) as the number of 
carbon atoms increased in the carbon chain. It is also well known, as 
discussed in Thaddeus (1994, 1995) and Vuitton et al. (2007), that 
unsaturated poiyynes can polymerize into long carbon linear chains 
from which fullerenes can form (Hunter, 1994; Thaddeus, 1994, 

1995) . As discussed in Vuitton et al. (2007) with photodissociation 
making the polyyne radical C 2 „H 2 +hv->-C 2n H+H which can then 
react with another polyyne C 2f ,H+C 2 mH 2 ->C 2 („ tm )H 2 +H to make a 
longer carbon chain, so very long carbon chains with very small 
amounts of H can form in Titan's atmosphere. The above reactions 
have been confirmed both experimentally and theoretically (Smith 
and Rowe, 2000; Smith et al., 2006; Berteloite et al., 2008). 
Laboratory measurements have also shown fullerenes can form 
large clusters (C 60 ) n with index n as high as n~146 (Martin et al., 

1996) . 

This does not mean to say that PAHs cannot also form in Titan’s 
atmosphere since the aromatic benzene C 6 H 6 and its radical 
phenyl C 6 H S are present in Titan’s atmosphere (Vuitton et al., 
2009a) and similar polymerization reaction paths may occur to 
make PAHs (Waite et al., 2007; Vuitton et al„ 2007). As we have 
stated previously both processes could very well be occurring at 
the same time and heavier aerosols could be mixtures or 
conglomerations of fullerenes and PAHs (see comments by 
Thaddeus, 1994, 1995 on large aerosol formation), 

2.8. Comment on UVIS observations of aerosols 

With regard to the UVIS observations by Liang et al. (2007) 
there could be a charge neutrality problem. Using Cassini UVIS 
stellar occupation data. Liang et al. (2007) observed the haze layer 
down to 400 km altitude. Assuming a scattering radius r A ~ 12.5 
ran for aerosols, they estimated a ''thoiin" particle density ~10 4 


part/cm 3 at 950 km and ~10 t> part/cm 3 at 400 km altitude. Such 
particle sizes could be in agreement with the Coates et al. (2007) 
heavy negative ion results if the particles are flat PAHs with radii 
r sce d<6nm. Since Liang et al.'s (2007) results are not very 
sensitive to particle size, UVIS could also be seeing ~6nm size 
or even smaller particles. The large number of "thoiin" particles 
observed by UVIS at 950 km ~10 4 part/cm 3 is much greater chan 
that observed by Coates et al. (2007) for negative ions and Waite 
et al. (2007) for heavy positive ions. 

if a significant fraction of the haze particles observed by UVIS 
are charged, then we have a conflict and charge neutrality would 
become a serious problem when the Langmuir Probe data is also 
considered for which N<, c <5000el/cm 3 (Wahlund, private com- 
munication, 2009). If Coates et al. (2007) are really seeing 
fullerenes then they may be too small (r A <12A) for UVIS to 
detect and Liang et al. (2007) are seeing particles too big to be 
observed by CAPS-ELS. But this would not resolve the charge 
neutrality problem unless all the haze particles observed by UVIS 
are not charged. For lower mass negative ions, as modeled by 
Vuitton et al. (2009b) and also found in ISM the ratio C 2 „H "7 
C 2 „H<?1, this ratio tends to increase with increase in the number 
of carbon atoms per carbon chain. Liang et al.’s (2007) results do 
show that the haze particles have a similar scale height to 
acetylene, which one might expect if the acetylene is the source of 
the “thoiin” particles. Acetylene can polymerize and make 
fullerenes while benzene is probably needed to make PAHs. This 
would favor more fullerenes relative to PAHs. 

2.9. Entrapped oxygen in aerosols 

Using the influx rate of keV of 1.1 x 10 6 O/cm 2 /s from Hartle 
et al. (2006b), a seed cross-section CT S eed~-4.5 x 10~ I4 cm 2 and 
N sre d~200 ions/cm 3 an implantation rate ~10~ 5 O implanted 
seed/cm 3 /s. The atmospheric volume needed to make one of 
the aerosols observed by Tomasko et al. (2005) is 
AV jonosph{ . re ~1.9 x 10 fi cm 3 , so there are ~19 O atoms/s implanted 
within each aerosol grain. If the seed particle fail rate is ~3 m/s 
and one uses a column density ~1001cm for 0 + penetration, we 
estimate ~6.3 x 10 5 O/aerosol. One can then estimate ~2 x 10 9 
fullerenes per aerosol, so that there are ~0.05% of the fullerenes 
with oxygen atoms. If one uses the medium mass flux limit by 
Borucki and Whitten (2008a, b) <P~4 x 10~ 13 kg/m 2 /s for the lower 
atmosphere one gets ~6.8% of the aerosols being in the form of 
fullerenes. But, these are lower limits since the heavy negative ion 
densities and masses are expected to increase below 950 km. 

The estimated mass flux of free oxygen to the surface is then 
0.1% of that incident onto the upper atmosphere from the 
magnetosphere, the bulk of the oxygen going to formation of CO 
and COj. 

Based on the above values for aerosol mass flux to the surface, 
the aerosol mass fraction from fullerenes, and fractional upper 
atmospheric oxygenation of the fuiierenes, we estimate ari oxygen 
mass flux to the surface of 3 x IQ -3 kg/m 2 per 300 Myr. This 
amounts to 7 x 10~ s mass fraction of the one-meter aerosol layer 
accumulated in 300 Myr. In the same accumulation time, abiotic 
atmospheric oxygen production condensed to the surface on 
Earth, mostly from UV photolysis of atmospheric water vapor is 
130 kg/m 2 (Chyba and Hand, 2001). For comparison, the mass 
fraction of photosynthetic oxygen in terrestrial sea water is 10" s 
(Schlesinger, 1997). 

2.10. Heavy ions and haze generation with oxygen model 

Fig. 4 summarizes the "heavy ions and haze generation with 

oxygen model" as presented in this paper. The top layer signifies 
Saturn’s magnetosphere with the exobase forming its lower 
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Heavy Ions and Haze Generation with Oxygen Model 
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Fig. 4 . Summary slide shows our “heavy ion and haze generation with oxygen 
model". See text for details. 


boundary. The ion composition of water group ions, methane 
group ions with energy l-10keV and more energetic H + , O * ions 
(10-100 keV). Below the exobase lies the topside ionosphere, with 
lower Ionosphere between 800 and 1400 km, where the keV 
oxygen ions can become trapped with fullerenes growing along 
with the PAH population. Cravens et al. (2008) show that the keV 
0 + ions deposit a significant fraction of their energy below 950 km 
altitude. 

We also show cosmic rays penetrating to the lower atmosphere 
and meteoritic dust down to low'er ionosphere. Below 950 km we 
assume further grain growth in size which then continues down 
to ~400 km. Below 200 km further growth presumably occurs at 
depths of maximum energy deposition from ionization by 
primary and secondary cosmic ray shower particles, in addition 
to grain growth, the cosmic rays can provide radiogenic energy to 
the aerosols which can then free the oxygen within the aerosols 
and possibly drive pre-biotic chemistry within the grains. Even 
after precipitation to the surface, the precipitating aerosol grains 
would continue to be processed chemically by irradiation from 
cosmic ray showers, also producing oxygen in surface or buried 
water ice to many meters in depth. Although height integrated 
cosmic ray ion production rates above 50 km maybe 100 times 
greater than that below 50 km (see Capone et al„ 1983; Molina- 
Cuberos et al., 1999), the time scales for descending aerosols 
~100-1000yr are much shorter than the surface irradiation times 
scales ~30-300 Mys (i.e., time it takes to accumulate ~lm of 
aerosols to surface). The descending aerosols may have a larger 
surface area for exposure than that on the surface. Therefore, 
cosmic rays will be effective in driving the aerosol chemistry at ail 
phases of their time history from formation at the highest 
altitudes -950 km down to the surface. 

Borucki et al. (2006) presented evidence for ion clusters such 
as CH|CH 4 , C7H7, C4H7 and H 4 C 7 N + to form below 400 km due to 
the ionization of nitrogen and methane by cosmic rays. They show 
that PAHs, which can be electrophilic (see Bakes et al., 2002), can 
form negative ions and significantly reduce electron densities, 
between 150 and 350 km, below that one would otherwise expect, 
in addition to PAHs fullerenes are also highly electrophilic as 
discussed previously and will similarly form negative ions. They 
also show that three-body electron-ion recombination collision 
will further reduce electron densities below expectations for 
heights <60 km (Huygens HAS1 conductivity probe measure- 
ments require low electron densities as discussed in Borucki and 
Whitten, 2008a,b). Borucki et al. (2006) also suggest that one 
must consider very small aerosols called embryos which can 
easily be ionized and can only retain charge states Q~l. These 
particles are similar in size to that observed above 950 km as 
heavy negative ions by Coates et al. (2007). Aerosol charging can 
contribute to atmospheric conductivities which can then be 
measured by the Huygens HAS! instrument (see Fulchignoni 


et al., 2005). The haze particles then fall through the troposphere 
to the surface. At the surface, due to accumulation of aerosol 
particles we show a mixture of "tholins", oxygen, nitrogen and 
amino acids with the latter a speculation. 

Below the surface a liquid ocean has recently been detected 
(Lorenz et al., 2008; Stiles et al., 2009) where interior models 
(Tobie et al., 2005) require ammonia as an anti-freeze since model 
temperatures are too low to allow pure liquid water to occur 
below the surface. It has been suggested that methanogenic life 
with this liquid ocean could be the source of Titan's methane 
which has a lifetime less than lOOMyrs (Fortes, 2000; McKay and 
Smith, 2005). This liquid water and heat, via the cryovoicanism 
that may have been observed (see Lopes et ai., 2007; Neish et al., 
2006), can then allow exobiological processes occur, especially if 
free oxygen is present. The observation of lakes by Cassini radar 
(Stofan et al., 2007), presents another tantalizing possibility for 
exobiological processes occurring at Titan. Recently, the Cassini 
VIMS experiment by Brown et al. (2008) has shown that these 
lakes contain ethane and are thus probably a mixture of organic 
molecules (see Raulin, 2008). If the aerosols discussed here do 
contain oxygen and land on these organic lakes a unique 
environment with significant astrobiological potential is created. 

How important could the magnetospheric oxygen be to the 
astrobiological chemistry of the Titan surface and liquid subsur- 
face environments? Over 300 Myr we estimate an accumulated 
fractional mass column density of 5-6 x 10 “ 7 of oxygen in a one- 
meter aerosol layer. This ratio would presumably remain 
representative of the exogenic material over time as this layer is 
mixed deeper into the regolith by geologic and impact processes. 
For comparison, the mass ratio of dissolved oxygen to water is 
10 5 in terrestrial sea water. This is also near the upper limits of 
radiolytic surface production from jovian magnetospheric inter- 
action with the surface of Europa for downward conveyance to its 
putative subsurface ocean (Cooper et al., 2001). By comparison, 
the exogenic oxygen at Titan is not significant at the macroscopic 
level for the bulk chemistry of the surface as compared to Earth 
and possibly Europa. 

On the microscopic level the astrobiological significance might 
be much higher if purely organic aerosol grains acquired ice shells 
at the surface, e.g. through water frost deposition from Enceladus- 
like cryovolcanic activity. The free oxygen abundance would have 
sources both from magnetospheric oxygen injection into the 
fullerenes and from direct radiolytic production in the water ice in 
all atmospheric and surface phases. As noted above, nitrogen and 
methane ions could also be implanted within the grains plus their 
addition to the aerosols as they fall through Titan’s atmosphere. 
Each Titan surface grain could be a microcosm of relatively 
abundant organics, water ice, and free oxygen trapped in 
fullerenes and the water ice. Heat from impacts and volcanic 
activity could then initiate the pre-biotic chemistry of life in 
microscopic versions of Darwin’s warm ponds on Titan. 
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